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We consider a model of active Brownian agents interacting via a harmonic attractive potential 
in a two-dimensional system in the presence of noise. By numerical simulations, we show that this 
model possesses a noise-induced transition characterized by the breakdown of translational motion 
and the onset of swarm rotation as the noise intensity is increased. Statistical properties of swarm 
dynamics in the weak noise limit are further analytically investigated. 



I. INTRODUCTION 



In different natural and social systems, agents form 
groups characterized by cooperative coherent motion p| . 
Such collective swarm motions have been observed and 
investigated in bacterial populations 12, I S H, IS B, in 
slime molds 0, IS IS IH3 1 f° r the ants lilt Irj and fish 
[l3L ItH Il5| , in the motion of pedestrians |16| and for the 
car traffic 17] . To describe these phenomena, various 
models of collective coherent motion in populations of 
self-propelled or self-driven particles have been proposed 
(see ^S El)- Some of them are formulated in terms 
of interacting automata located on a lattice or having 
continuous coordinates [IS El US ■ A different class 
of models is based on dynamical equations for individ- 
ual self-propelled particles, including velocity-dependent 
negative friction terms to account for their active mo- 
tion and assuming that interactions between the parti- 
cles are described by some binary potentials. This latter 
approach has been used to phenomenolqgically charac- 
terize motion of biological organisms 124). |25l. individ- 
ual blood cells |M E3 and humans ll6L l28i . and to de- 



scribe the behavior of physical particles with energy de- 
pots [2S, H3, EIH HS| • Effective dynamical equations with 
velocity-dependent friction could be approximately de- 
rived for floating particles which propel themselves by 
releasing surfactant into the medium |33| . A contin- 
uum approximation to discrete particle models, hydro- 
dynamical models of active fluids have been proposed 

ESIIIIISIMIIIIII. 

The common property of all swarm models is that they 
show the emergence of coherent collective flows starting 
from a disordered state with random velocity directions 
of individual particles (agents) . Such kinetic transitions 
have been extensively investigated for the automata sys- 
tems I2S l39j and in the framework of hydrodynamics 
[34L l37j|. The ordered states of swarms can represent 
simple translational motion or be characterized by vor- 
tex flows. Both spatially distributed and localized swarm 
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states are possible. In a distributed state, the population 
fills the entire available medium. In contrast to this, the 
population forms compact spatial groups in a localized 
state. An interesting example of a localized swarm state 
is a rotating flock of finite extent, seen in the simula- 
tions of a discrete model of self-propelling particles and 
described analytically within a continuum active fluid ap- 
proximation [2SI (see also |T^|). 

Localized states of swarms may undergo transitions 
leading to new dynamical regimes, when the system pa- 
rameters or the noise intensity are gradually changed. 
In the previous publication |4C|. a noise-induced transi- 
tion from the localized state with translational motion 
to a state with incoherent oscillations without transla- 
tional motion was investigated numerically and analyt- 
ically for a one-dimensional system of interacting self- 
propelled particles. In the present article, we extend in- 
vestigations of this system to two spatial dimensions. We 
study here a population of self-propelled particles inter- 
acting via a parabolic interaction potential corresponding 
to linear attracting forces between the pairs. In absence 
of noise, this dynamical system has two kinds of attrac- 
tors, corresponding, respectively, to a compact traveling 
state of the entire population and to a state where it ro- 
tates as a vortex without any global translational motion. 
The aim of our study is to investigate effects of noise on 
translational swarm motion. We find that the system 
is highly sensitive to stochastic forces. When noise is 
present, the traveling swarm is a cloud of particles char- 
acterized by different dispersions in the directions parallel 
and transverse to the direction of translational motion. 
Our numerical simulations confirmed by an approximate 
analytical study show that the mean-square transverse 
dispersion of a swarm is proportional to the square root 
of the noise intensity, whereas its dispersion along the di- 
rection of motion depends linearly on the noise intensity. 
Therefore, for weak noises the swarm looks like a pancake 
oriented orthogonally to the motion direction. When the 
noise is increased, the swarm gradually acquires a more 
symmetric shape. For strong noises, we find that the 
translational motion of a swarm becomes suddenly im- 
possible and is abruptly replaced by a rotational regime 
with a vortex flow. The detailed formulation of the model 
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is presented in the next section. In Sect. Hill we describe 
the results of numerical simulations. The statistical prop- 
erties of a traveling swarm in the weak noise limit are ap- 
proximately explained by an analytical theory which is 
constructed in Sect. II Vl The paper ends with conclusions 
and discussion of obtained results. 



II. THE MODEL 

We consider a swarm formed by N identical self- 
propelled particles of unit mass interacting via an attrac- 
tive parabolic pair potential. The dynamics of the system 
is given by the following set of evolution equations: 



(la) 
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The forces Fj depend on the particle velocity and are 
introduced to take into account active motion of particles. 
We choose them in the form 



Fi = (1 - v? 



(2) 



so that, in absence of noise and interactions, the particle 
acquires the unit velocity v = 1. Additionally, the parti- 
cles are subject to stochastic white forces & of strength 
D which are independent for different particles and are 
characterized by the correlation functions 

«<(*)> =o, m)^(t , ))=2DS(t-t% j . 

This model has previously been introduced in jljj. It 
is phenomenological, but rather general because it can be 
viewed as a normal form for a population of particles near 
a supercritical bifurcation corresponding to spontaneous 
onset of active motion (see 19]). In this model, attrac- 
tive interactions between particles have an infinite range 
and grow linearly with the separation between them. Be- 
cause we shall consider only spatially localized states of 
the population, our results will hold, however, also for 
the situations when interactions are characterized by a 
finite range, but it is much larger than the mean swarm 
diameter. It should be noted that, in a different context, 
the model Q has been considered already by Rayleigh 
0. 

The study of the one-dimensional version of the 
model Q has shown that, as the noise intensity D is in- 
creased, spontaneous breakdown of translational swarm 
motion takes place here |4jj. Some statistical proper- 
ties of translational swarm motion in the two-dimensional 
model (1) (with a slightly different choice of the forces 
have subsequently been discussed |43- For the case of 
two interacting particles (N = 2), the rotational modes 
were described in |32j . where simulations for small ro- 
tating clusters consisting of 20 particles have also been 
reported. The aim of the present work is to perform sys- 
tematic, numerical and analytical, investigations of the 
behavior described by this two-dimensional model. 



III. NUMERICAL SIMULATIONS 

Numerical integration of equations Q was performed 
using the Euler scheme with the constant time step of 
0.001. In all simulations, the total number of particles 
was fixed to N = 300 and the coefficient a specifying 
the strength of interactions between the particles was set 
to a = 100. To produce a traveling localized state of 
the swarm, special initial conditions were used. At time 
t = 0, all particles had identical positions and velocities 
and the noise was switched only a little later, at time 
t = 30. 

Several statistical characteristics of the swarm were 
monitored during simulations. The center of mass R of 
the swarm and its mean velocity V at time t were defined 
as R(t) = {I /N) Ei r ( (t) and V(t) = (1/N) Vi (i), re- 
spectively. Because the cloud of traveling particles in the 
presence of noise was anisotropic, we also determined its 
instantaneous mean-square dispersions in the directions 
parallel (Sn) and orthogonal (S±) to the direction of its 
instantaneous mean velocity V. They were defined as 

1 N 

S\l(t) = ^™yE[(r*W-RW)'V(t)] 2 (3a) 

1 N 

Sx(t) = ^^£[(r t (i)-R(i))xV(t)] 2 .(3b) 

^ ' i— 1 

Additionally, angular momenta of all particles with re- 
spect to the mass center were determined as 



Li(t) = (r t -(t) - R(t)) x ( Vi (t) - V(i)) • 



(4) 



Figure 
tude V = 



la 



shows the time dependences of the magni- 
|V| of the mean swarm velocity for two simu- 
lations with different noises. When the noise is relatively 
weak (D = 0.067), its introduction leads to some fluctua- 
tions in the instantaneous swarm velocity and a decrease 
of its average level. If a stronger noise (D = 0.07) is ap- 
plied, the swarm velocity first behaves as for the weaker 
noise, but then abruptly drops down to a value close to 
zero. This sudden transition corresponds to the break- 
down of translational motion of the swarm. In Fig. 1(b) 
we have plotted the time-averaged swarm velocity (V) as 
the function of the noise intensity D. The average veloc- 
ity gradually decreases with noise, until the breakdown 
occurs at 0.067 < D < 0.070. 

In the state with translational motion at relatively 
weak noises, the direction of the swarm motion does not 
remain constant with time. The swarm travels along a 
complex trajectory, part of which is shown in Fig. [2 (such 
trajectories should correspond to the Brownian motion of 
the entire swarm). In the inset in this figure, we display 
the distribution of particles in the swarm at some time 
moment. It can be noticed that the cloud of particles is 
significantly squeezed along the direction of swarm mo- 
tion. 

Figure[3]shows the computed average longitudinal (Sh ) 
and transverse (S±) dispersions as functions of the noise 
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FIG. 1: (a) Time evolution of the mean velocity of a swarm 
before (solid line, D — 0.067) and after (dashed line, D = 
0.070) the stochastic breakdown of the translational mode, 
(b) Mean velocity of the swarm with increasing noise strength. 
At a critical noise strength a sharp transition in the behavior 
of the swarm occurs (see also Fig. 



intensity D. For weak noises, S± 3> Su so that the swarm 
is strongly squeezed. As noise increases, the shape of 
the swarm becomes more symmetric and the transver- 
sal dispersion approaches the dispersion along the direc- 
tion of translational motion. Finally, after the break- 
down of translational motion has taken place for a suffi- 
ciently strong noise, the swarm becomes statistically cir- 
cular {S± = 

The sequence of snapshots in Fig. 0] displays temporal 
evolution of the swarm when the noise intensity exceeds 
the breakdown threshold. Initially, the swarm is travel- 
ing and its shape is similar to that characteristic for the 



FIG. 2: Motion of the center of mass of the swarm (300 par- 
ticles) within a certain time window for noise below a critical 
one. For t = 1140 the corresponding snapshot of the swarm 
is shown. The red arrow shows here th mean swarm velocity; 
the noise intensity is D = 0.067 
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FIG. 3: Behavior of the longitudinal (solid line [+]) and 
transversal (dashed line [x]) dispersions of the swarm with 
increasing noise. 



weaker noises (cf. Fig.EJ. However, in the course of time 
the swarm slowly acquires a ring shape, with particles 
rotating around its center. This rotating ring structure 
corresponds to a state where translational motion of the 
entire swarm is already absent. 

A different visualization of the process, accompanying 
the breakdown of translational motion and the transi- 
tion to a rotating swarm, is chosen in Fig. EJ Here we 
show the trajectory of motion of the center of mass of the 
swarm (solid line) together with the trajectory of motion 
of one of its particles (dashed line). We see that, in a 
traveling swarm, the particles perform irregular oscilla- 
tions around its instantaneous mass center. When the 
translational motion is terminated and the rotating ring 
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FIG. 5: The trajectories of the center of mass (solid line) and 
of a single particle are shown (dashed line) above the critical 
noise strength (for D = 0.070). 

are present in the considered model. It can be expected 
that, if short range repulsive interactions arc additionally 
introduced, the breakdown of the rotational symmetry in 
the ring would occur and one of the rotation directions 
would be selected [82^ . 




FIG. 4: Several subsequent snapshots of a swarm with 300 
particles during the transition from translational motion to 
the rotational mode; the noise intensity is D — 0.070. 



is formed, such oscillations become transformed into ro- 
tations around the ring center. 

To provide statistical description of particle motions 
in the traveling and rotating states of the swarm, angu- 
lar momentum distributions P(L) have been computed. 
For the state with translational motion, the distribution 
has a single central peak at L — (Fig. |6(a)| ). In con- 
trast to this, in the rotational state the distribution has 
two symmetrically placed peaks corresponding to a cer- 
tain non- vanishing momentum (Fig. 6(b) I. Note that the 
particles inside the ring are rotating both in the clock- 
wise and counter-clockwise directions, and the numbers 
of particles rotating in each direction are approximately 
equal. Thus, the swarm does not rotate as a whole and 
its total angular momentum remains zero on the aver- 
age. This behavior is the consequence of the fact that 
only long-range attractive interactions between particles 



(a) D = 0.067 



(b) D = 0.070 



FIG. 6: Distribution of angular momenta of the particles (a) 
in the traveling (D = 0.067) and (b) rotating (D = 0.070) 



IV. THE WEAK NOISE LIMIT 

Our numerical simulations have shown that, for weak 
noises, the swarm is strongly squeezed in the direction 
along its center-of-mass motion and its longitudinal and 
transverse dispersions are strongly different. Below we 
derive approximate analytical expressions for Su and 
Sjin the limit of the small noise intensity D — ► 0. 

First, we note that in this limit the motion of the center 
of mass of the swarm will remain approximately linear 
within very long times or, in other words, the swarm 
velocity V remains approximately constant on the short 
time scales characteristic for the motions of individual 
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particles inside the traveling swarm. We introduce the 
coordinate system in such a way that its x-axis is parallel 
to the direction of the swarm motion and its y-axis is 
orthogonal to it. The coordinates Xi and yi of all particles 
forming the swarm can be written as Xi = X + Sxi and 
i/i = Y + Syi where X and Y are the coordinates of the 
swarm center R. By our choice, we have Y = 0, so that 
Vi = Syi. 

To derive the evolution equation for X, we notice that 



*(*) = ^$>(t) = <*i(*)> 



(5) 



Summing up the evolution equations for all Xi, we ap- 
proximately obtain 



X- (l-X 2 )X-3(5x t )X-(y 2 )X = (6) 

where we have neglected the terms with the higher pow- 
ers of velocity fluctuations Sxi and yi. 

In the statistical steady state, X = V = const and © 
is reduced to the equation 



(1-V 9 )- 3 «)-<»?) 



V = 



(7) 



determining the velocity of swarm motion in the presence 
of noise. Its solution for the traveling swarm (V ^ 0) is 



V 2 



1 



■(&!)-<#> 



(8) 



The evolution equation for 5xi can be obtained by sub- 
tracting equation JBJ from the equation for the variable 
Xi in the model Keeping only the leading terms, 

linear in deviations from the mass center, we get 



Sxi + 2Sxi + aSxi = £ x (t) 



(9) 



This is an evolution equation for a damped harmonical 
oscillator. Note that fluctuations of Xi(t) are not coupled 
to the transverse component yi(t). 

In a similar way, the evolution equation for the trans- 
verse deviations yi{t) can be obtained, 

Hi - (1 - V 2 ) Vl + 2VSx iyi + yf + ay, = £»(t) (10) 

In this equation, we have retained nonlinear terms. This 
is done because such terms are essential for the damping 
of oscillations of the transverse component. 

Indeed, if such terms are neglected, we would have 



y l -(l-V 2 )y l +ay l =^(t) 



(11) 



Because, as follows from equation |JS}, we have V 2 < 
1, oscillations in yi would then exponentially grow with 
time. Thus, nonlinear terms play a principal role for 
transverse fluctuations and cannot be neglected even in 
the weak noise limit, in contrast to the respective terms 
for the longitudinal fluctuations 5xi. 



As shall be verified at the end of our derivation, the 
condition (yf) ^> (Sx^ holds in the weak noise limit. 
Therefore, the swarm velocity is mostly influenced by 
the transverse fluctuations and we have approximately 
V 2 = 1 - (y 2 ). Substituting this into l(TD)l. we get 



(12) 



The stochastic differential equation (|12|) does not in- 
clude longitudinal fluctuations 5xi and, furthermore, 
fluctuations for different particles i are not coupled here. 
For subsequent analysis of this differential equation, we 
drop the indices and write as 



v- ((f) - f) y + ay = m 



(13) 



with (t)) = 2DS(t - t'). Here we have assumed 

that statistical averaging is equivalent to averaging over 
the ensemble. 

The approximate solution for the probability distribu- 
tion of variable y in the statistical stationary state can 
be derived for this equation assuming that the parameter 
a, determining the oscillation frequency, is large (a» 1). 
We introduce slowly varying amplitudes 



y(t)=r ] (t)e i " t +r)*(t)e- 



(14) 



where u> — yfa ^S> 1. Substituting this into Ijl3(l and keep- 
ing only the resonant terms of the highest order inw, we 
obtain a stochastic evolution equation for the transverse 
complex oscillation amplitudes 



n 



\v\ 



„--u, 2 M 2 „ + C(t) (15) 



where the complex-valued white noise has correlation 
functions 



(C(t))=Q, (CWC(O> = 0, (C(t)C(t')) 



D 
2^2 



S(t-t') 



This stochastic Langevin equation corresponds to the 
following Fokker-Planck equation for the probability den- 
sity P = P(r),r)*,t): 



OP 
dt 



d_ 
dr] 
_d_ 

cV 
1 



"[{\ri\ 2 )'n-~\4 



r)*P 



(16) 



d 2 P 

2oj 2 dr/drj* 

The stationary solution P of the Fokker-Planck equation 
reads 



P = — exp 



-^(-4(M 2 )M 2 + 3M' 



(17) 



where the normalization constant Z is given by 

,4 



Z = I exp 



--(-4 
D 



V\ 



I "'/I 



3H 



tfi). (18) 
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When the probability distribution is known, the second 
statistical moment can be calculated as 



<H 2 ) 



\v 



2 P(r,,r,*)d 2 r ] 



(19) 



Because the stationary probability distribution (jT7|) de- 
pends on (|r/| 2 ) , this is an equation which should be 
solved to determine this statistical moment. Let us 



substitute 77 



pe 



^ and uj 2 = 



Then 



= p and 



d 2 r] = pdpdcj). Therefore, equation 119fl takes the form 



r OO Q 

Jo P ex P 


-^(-4<P 2 )P 2 




dp 


r °°pexp 


-g(-4<p»)pM 


-3p 4 )]dp 



(p 2 ) 



Introducing the variable u 
is transformed to 



(20) 



p(H 2 ) 1/2 , equation fZ® 



1 



/ °° u 3 exp [-v (-4u 2 + 3u 4 )]du 
/ °° wexp [— v (— 4« 2 + 3u 4 )]du 



where 



(21) 



(22) 



Numerical solution of equation l|21|) yields v ~ 0.22. 
When ^ is known, equation l|22|) determines (|??| 2 ) as 



(\r ] \ 2 )=v 1 l 2 B 1 l 2 a- 1 . 



(23) 



Using the definition (|14|) of variable rj, we find that 5^ = 
(y 2 ) = 2(|r/| 2 ). Thus, we finally obtain the analytical 
estimate for the transverse dispersion of the swarm in 
the weak noise limit, 



Si 



£1/2 



(24) 



where the numerical coefficient is k = 2^ : / 2 ~ 0.94. 

The longitudinal dispersion S\\ = (<5x 2 ) is approxi- 
mately determined by the linear stochastic differential 
equation ©. A straightforward derivation for a ^ 1 
yields for this property the analytical estimate 



S\\ = 



D 
2a' 



(25) 



We see that in the limit D — > the longitudinal dis- 
persion is indeed much smaller than the transverse dis- 
persion of the traveling swarm, as assumed in the above 
derivation. 

Note that, for a ^> 1, statistical dispersions of trans- 
verse and longitudinal velocity fluctuations W± = (y 2 ) 

and W\\ = (^5 x 2 ^ are W± = aS±_ and W\\ = aS\\. There- 
fore, in the considered weak noise limit they are given 
by Wj_ = K-D 1 / 2 and Wu = D/2. Comparing these ana- 
lytical estimates with the simulations in Fig. [7| we find 
that they agree well with the respective numerical data. 
For the transverse dispersion, the agreement is found in a 
wide interval of noise intensities. The analytical expres- 
sion for the longitudinal dispersion holds, on the other 
hand, only for very weak noises. 
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FIG. 7: Dispersions of swarm velocities in (a) transversal and 
(b) longitudinal directions as functions of the noise intensity. 
The symbols show the simulation data. The dotted lines are 
the theoretically predicted power law dependences. 



V. CONCLUSIONS 

We have studied statistical properties of localized 
swarms with long-range attractive interactions in two- 
dimensional media. Our numerical simulations show that 
the swarm is highly sensitive to the action of noise. Even 
very weak noises lead to strong dispersion of the swarm 
along the direction orthogonal to the direction of its 
translational motion. The approximate analytical the- 
ory predicts that the transverse dispersion of aswarm in- 
creases as vT> with the noise intensity D, whereas its 
longitudinal dispersion depends linearly on D in the limit 
D — ► and remains therefore much smaller in this limit. 
Hence, for weak noises the traveling swarm is strongly 
squeezed along the direction of its mass motion. This 
analytical result is confirmed by numerical simulations. 

Increasing the noise intensity D, we find that transla- 
tional motion breaks down when a certain critical inten- 
sity is reached. After the breakdown, the translational 
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motion is stopped and instead the swarm goes into a rota- 
tional mode where the center of mass of the swarm shows 
only weak random motion. This behavior resembles the 
breakdown of translational motion which was previously 
seen for the one-dimensional system |4fj. In contrast to 
the one-dimensional case, we could not however analyti- 
cally treat this transition, because of the strong fluctua- 
tions in the transverse direction. 

Though our results are obtained in the model with har- 
monic attractive interactions, they are also be applicable 
for the models with finite-range attractive interactions 
between the particles, provided that the size of a localized 
swarm (i.e., the statistical dispersion of the coordinates 
of its particles with respect to the mass center) is much 
smaller than the interaction radius. In this situation, a 
harmonic approximation of the interaction potential can 



be applied. Moreover, similar effects can be expected for 
swarms in three-dimensional media. It would be inter- 
esting to see whether the discussed behavior is indeed 
characteristic for real biological swarms with long-range 
interactions between individual organisms, such as bird 
flocks or fish schools. 



Acknowledgments 

We are grateful to L. Schimansky-Geier and D. Zanette 
for useful discussions. This study has been performed 
in the framework of the Collaborative Research Cen- 
ter "Complex Nonlinear Processes" of the Deutsche 
Forschungsgemeinschaft (DFG-SFB 555). 



E. Ben-Jacob, I. Cohen, and H. Levine, Adv. Phys. 49, [23] 
395 (2000). 

R. M. Harshey and T. Matsuyama, Proc. Nat. Acad. Sci. [24 
USA 91, 8631 (1994). [25 
E. O. Budrene and H. C. Berg, Nature 376, 49 (1995). [26 
E. Ben- Jacob, I. Cohen, A. Czirok, T. Vicsek, and 
D. Gutnick, Physica A 238, 181 (1997). [27 
M. P. Brenner, L. S. Levitov, and E. O. Budrene, Bio- 
phys. J. 74, 1677 (1998). [28 
A. Czirok, M. Matsushita, and T. Vicsek, Phys. Rev. E [29 
63, 031915 (2001). 
H. Levine and W. Reynolds, Phys. Rev. Lett. 66, 2400 [30 
(1991). 

D. A. Kessler and H. Levine, Phys. Rev. E 48, 4801 [31 
(1993). 

S. Nagano, Phys. Rev. Lett. 80, 4826 (1998). [32 
W.-J. Rappel, A. Nicol, A. Sarkissian, H. Levine, and 
W. F. Loomis, Phys. Rev. Lett. 83, 1247 (1999). [33 

E. Bonabeau, G. Theraulaz, V. Fourcassie, and J.-L. 
Deneubourg, Phys. Rev. E 57, 4568 (1998). 
A.-C. Mailleux, J.-L. Deneubourg, and C. Detrain, Anim. [34 
Behav. 59, 1061 (2000). [35 
G. Flierl, D. Griinbaum, S. A. Levin, and D. Olson, J. [36 
theor. Biol. 196, 397 (1999). 

J. K. Parrish and L. Edelstein-Keshet, Science 284, 99 [37 
(1999). [38 
J. K. Parrish, S. V. Viscido, and D. Griinbaum, Biol. [39 
Bull. 202, 296 (2002). 

D. Helbing and P. Molnar, Phys. Rev. E 51, 4282 (1995). [40 

D. Helbing, Rev. Mod. Phys. 73, 1067 (2001). 

A. Okubo and S. A. Levin, Diffusion and Ecological Prob- [41 
lems: Modern Perspectives (Springer, New York, 2001), 
2nd ed. [42 
A. S. Mikhailov and V. Calenbuhr, From Cells to Soci- 
eties (Springer, Berlin, 2002). [43 
T. Vicsek, A. Czirok, E. Ben- Jacob, I. Cohen, and 
O. Shochet, Phys. Rev. Lett. 75, 1226 (1995). 

E. V. Albano, Phys. Rev. Lett. 77, 2129 (1996). 
A. Czirok and T. Vicsek, Physica A 281, 17 (2000). 



H. Levine, W.-J. Rappel, and I. Cohen, Phys. Rev. E 63, 
017101 (2001). 

H.-S. Niwa, J. theor. Biol. 171, 123 (1994). 

H.-S. Niwa, J. theor. Biol. 181, 47 (1996). 

M. Schienbein and H. Gruler, Bull. Math. Biol. 55, 585 

(1993). 

M. Schienbein and H. Gruler, Phys. Rev. E 52, 4183 
(1995). 

D. Helbing, Verkehrsdynamik (Springer, Berlin, 1997). 
F. Schweitzer, W. Ebeling, and B. Tilch, Phys. Rev. Lett. 
80, 5044 (1998). 

W. Ebeling, F. Schweitzer, and B. Tilch, BioSystems 49, 
17 (1999). 

U. Erdmann, W. Ebeling, F. Schweitzer, and 
L. Schimansky-Geier, Eur. Phys. J. B 15, 105 (2000). 
U. Erdmann, W. Ebeling, and V. S. Anishchenko, Phys. 
Rev. E 65, 061106 (2002). 

A. S. Mikhailov and D. Meinkohn, in Stochastic Dy- 
namics, edited by L. Schimansky-Geier and T. Poschel 
(Springer, Berlin, 1997), pp. 334-345. 
J. Toner and Y. Tu, Phys. Rev. Lett. 75, 4326 (1995). 
Z. Csahok and A. Czirok, Physica A 243, 304 (1997). 

U. Erdmann, Interjournal, 114 (1997), URL 

http: //www. interjournal . org/manuscript_abstract . php? 10536 

J. Toner and Y. Tu, Phys. Rev. E 58, 4828 (1998). 

U. Erdmann, Kollektive Bewegung (Logos, Berlin, 2004). 

A. Czirok, A.-L. Barabasi, and T. Vicsek, Phys. Rev. 

Lett. 82, 209 (1999). 

A. S. Mikhailov and D. Zanette, Phys. Rev. E 60, 4571 
(1999). 

J. W. S. Rayleigh, The Theory of Sound, vol. I (Dover, 
New York, 1945), 2nd ed. 

F. Schweitzer, W. Ebeling, and B. Tilch, Phys. Rev. E 
64, 021110 (2001). 

U. Erdmann, in Proceedings from the International Con- 
ference on Complex Systems on Unifying Themes in 
Complex Systems, edited by Y. Bar- Yam (Perseus Books, 
Cambridge, Massachusetts, 2000), pp. 153-161. 



